Flexoelectricity and piezoelectricity - reason for rich variety of phases in 

antiferroelectric liquid crystals 
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The free energy of antiferroelectric liquid crystal which takes into account polar order explicitly is 
presented. Steric, van der Waals, piezoelectric and flexoelectric interactions to the nearest layers 
and dipolar electrostatic interactions to the nearest and to the next nearest layers induce indirect tilt 
interactions with chiral and achiral properties, which extend to the third and to the fourth nearest 
layers. Chiral indirect interactions between tilts can be large and induce helicoidal modulations 
even in systems with negligible chiral van der Waals interactions. If indirect chiral interactions 
compete with chiral van der Waals interactions, the helix unwinding is possible. Although strength 
of microscopic interactions change monotonically with decreasing temperature, effective interlayer 
interactions change nonmonotonically and give rise to nonmonotouous change of modulation period 
through various phases. Increased enatiomeric excess i.e. increased chirality changes the phase 
sequence. 



PACS numbers: 61.30.Cz 

Chiral polar smectics are materials which form lay- 
ered phases. In these phases the average direction of 
elongated molecules is tilted for an angle with respect 
to the layer normal. Smectic layers are spontaneously 
polarized in the plane of the layer and perpendicular to 
the direction of the tilt |l|] . Among these polar smectics 
the most widely experimentally and theoretically studied 
materials are ferroelectric liquid crystals. Ten years ago 
at an attempt to synthesise ferroelectric liquid crystals 
with larger polarization, antiferroelectric liquid crystals 
were discovered [^J. Their synthesis stimulated exten- 
sive experimental research which resulted in a discovery 
of many new polar liquid crystalline phases. In addi- 
tion to the the ferroelectric Sm C* phase and the an- 
tiferroelectric Sm phase, the Sm C* phase and the 
Sm C* phase were found to exist in some materials while 
in others instead of the Sm C* phase, two phases - the 
Sm C* FI2 phase and the Sm C* FI1 phase are present. The 
phases are modulated along layer normal and the mod- 
ulation period of these phases was determined by reso- 
nant X-ray scattering || . The period was shown to vary 
nonmonotonically with decreasing temperature from a 
few layers in the Sm C* phase, to a few hundred lay- 
ers in the Sm C* phase, to approximately four layers in 
the Sm C* FI2 phase, to approximately three layers in the 
Sm C^ 71 phase and to approximately two layers in the 
Sm phase at low temperatures. 

Antiferroelectric liquid crystals can be theoretically 
described by continuous models Q-g that predict only 
two of experimentally observed phases, the ferroelectric 
Sm C* phase and the antiferroelectric Sm phase. Dis- 
crete models |7[|| take interlayer interactions between 
the tilt vectors explicitly into account. The model 
considers competing interactions between tilts in near- 



est layers and in next nearest layers and predicts various 
structures where the direction of the tilt varies uniformly 
from layer to layer. The model was later called the clock 
model |J. 

Within this model the observed nonmonotonous tem- 
perature dependence of the modulation period can be re- 
produced only with nonmonotonous temperature depen- 
dences of model parameters, which cannot be understood 
from microscopic interactions. Although the model can 
predict structures of the Sm C* FI2 and of the Sm C^ n for 
the specific set of clock model parameters, these phases 
cannot exist as stable phases in a broader temperature 
interval. For this some synclinic interactions would be 
needed between the third and between the fourth near- 
est layers ||. 

In this letter we present the model where the polar- 
ization of smectic layers is taken explicitly into account 
and is treated as a secondary order parameter to the tilt. 
In addition to the steric and van der Waals interactions 
to the nearest layers between the tilt vectors we consider 
also piezoelectric and flexoelectric couplings between the 
tilt and the polarization and electrostatic interactions be- 
tween the polarizations of the nearest layers and polar- 
izations with the next nearest layers |hJ. 

To describe the structure of phases, we introduce the 
tilt vector £j which gives the magnitude and the direction 
of the tilt in the j-th smectic layer and the polar order 
parameter r)j = {rjj X ,r]jy} which gives the magnitude 
and the direction of the in plane transverse polar order 
in the j-th layer. The spontaneous polarization in the j- 
th layer is proportional to the polar order tj: Pj = Pq r/j, 
where Pq is a layer polarization for a completely polarly 
ordered layer . Therefore polar order parameter r\ will 
be here shortly called •polarization. The free energy of the 
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smectic system with N layers can be expressed in terms 
of both order parameters as: 

N 1 1 1 

G = 12 a Go + 2 ai & ' + 2 fl ^ j X + 

3=1 

+ J & o »7 3 2 + ^ &i faj • Vj+l) + l b 2 (Vj ■ V3+2) + 
+ c p (Vj x + (£3+1 - €3-1) • (i) 

Here G = 1/2 o(T - T ) + 1/4 6 £4 is the tilt depen- 
dent free energy of an isolated layer, where a continuous 
transition to the tilted phase takes place at To. Divid- 
ing all terms by the parameter a, we express all model 
parameters as well as the free energy G in the units of 
temperature. Parameter a,\ corresponds to steric and van 
der Waals interactions between nearest layers. Steric in- 
teractions favor synclinic molecular alignment in neigh- 
boring layers and give a negative contribution to ax, while 
van der Waals intermolecular interactions favor anticlinic 
ordering in neighboring layers and give a positive contri- 
bution to a\. The parameter a\ is of the order of a few 
Kelvins @Jll]] . We believe that steric and van der Waals 
interactions are negligible further than to the nearest lay- 
ers. The parameter /1 gives the chiral part of van der 
Waals interactions between molecules |l2| in neighbor- 
ing layers and vanishes in racemic mixtures. We assume 
that /1 is smaller than a\ and it is a few tenths of a 
Kelvin. Polar vectors are associated with the layer polar- 
ization via magnitude and direction of molecular dipole 
moments. The parameter 60 is positive because polariza- 
tion is not a proper order parameter and both b\ and 62 
favor antiparallel dipolar order and are therefore positive 
f[o|| . Since they give intralayer and interlayer interac- 
tions which decrease for more distant layers we expect 
that 60 >> b\ >> 62. Transition temperatures from the 
Sm C* phase to the Sm phase are few degrees lower 
in nonpolar racemic mixtures than in pure samples |^] 
which suggests that 6 is also of the order of a few K. 
Polarization is induced in chiral polar smectics by the 
tilt and the piezoelectric coupling is given by c p . The 
parameter c p has also a chiral character and vanishes in 
racemic mixtures. The magnitude of the layer polariza- 
tion is only one tenth of the polarization of completely 
polarly ordered system [jl3| i.e. the magnitude of rjj is 
of the order 0.1 which means that also the value of c p is 
of the order of a Kelvin. The polarization coupling with 
the tilt in neighboring layers is given by /1 and is of the 
same order as c p |l~4| , |l5| ] . 

The polar part of the free energy Eq. (Q), can be writ- 
ten in the form 

Gp=vC£ + ±r,B V . (2) 

The tilt vector is 2iV-dimensional vector of the form 
€ = {fix, —£jx, ■■4Nx,£,iy,--£,Ny} and similar for the 



polarization r). The five-diagonal 2iV-dimensional ma- 
trix B_ gives intralayer and interlayer electrostatic in- 
teractions with elements Bjj — bo,Bjj±i = i b% and 
Bj,j±2 — \ ^2- The 27V-dimensional matrix C_ gives 
the coupling between tilts and polarizations. The off- 
diagonal elements Cjj±i = Cj+N,j+N±i = n H an d 
Cj,j+N = ~Cj+N,j = Cp- Minimization of Eq. (Q) with 
respect to polarization r\ gives 

r) = -B- 1 ■ C£, (3) 

where is the inverse matrix of the five-diagonal 

matrix B_. As 6 » b\ » 62 we keep in B _1 
only the terms up to the second order in 61/60 and up 
to the first order in 62/60. The elements of B^ 1 are 

^ = + i t^f) = (fe) 

BTj ±2 = £ (2 (^) 2 - (I)). Inserting the tilt dc- 

pendent polarization Eq. ([}]) into Eq. (^) we obtain the 
polar part of free energy as 

G P = -\tiCB- 1 Cti (4) 
and the free energy due to interlayer interactions is 

1 N ( 4 3 \ 

Ant = 2 12 [12 hk (& ■ + 12 ? k & x ^'+ fe) • 
j=l \k=l fe=l / 

(5) 

Parameters hk and fk which appear after polarization 
elimination, give effective interactions between tilts. Al- 
though direct van der Waals and steric interactions are 
significant only to the nearest layers and electrostatic in- 
teractions are significant up to the next nearest layers, 
effective interactions are significant up to the fourth near- 
est layers. 
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Achiral effective interactions between nearest layers a\ 
consist of direct steric and van der Waals interactions (ai) 
as well as of indirect interactions due to piezoelectrically 
and flexoelectrically induced polarization. Achiral inter- 
action between next nearest layers are only indirect and 
can be for systems with negligible flexoelectric interaction 
(fi ~ 0) either competing (a 2 > 0) or noncompeting 
with (5,2 < 0) [ p^| . The second case leads to the validity 
of the continuous bilayer models ||-[| , although it seems 
that the present experimental knowledge supports this 
possibility by a single experimental evidence fllq] . In sys- 



tems with large flexoelectric interactions (|c p 



5,2 is always positive and competes with interactions be- 
tween nearest layers 5\. Indirect effective interactions 
between third neighbors are given by 0,3. The term is 
always negative and favors synclinic tilt directions in in- 
teracting layers. Therefore this term in systems with 
significant flexoelectric interactions tends to stabilize the 
structures with three layer periodicity. Indirect effective 
interactions between fourth nearest layers 04 can be ei- 
ther negative or positive. 

The effective chiral coupling between neighboring lay- 
ers has two contributions: the van der Waals originated 
(/1) and the polarization part which appears due to the 
combination of piezo and flexoelectricity (c p n/bo). Even 
systems with negligible direct chiral interactions f± ~ 
are hellicoidally modulated. In all previous models, chi- 
ral interactions have always been treated as weak interac- 
tions which induce only slight perturbations in the struc- 
ture i.e. hclicoidal modulations with very long periods. 
In antiferroelectric liquid crystals this is not always the 
case. The ratio c p [ijb^ can have the value of a few 
Kelvins and is comparable to other achiral interlayer in- 
teractions. The competition between the van der Waals 
part of chiral interactions (/1) and indirect chiral interac- 
tions due to piezoelectric and flexoelectric effect (c p /V&o) 
can explain the helix unwinding without polarization re- 
versal as observed in the antiferroelectric phase of various 
systems ||[I6). 

The simplest solution, which minimizes the free energy 
Eq. (|J) for competing interlayer interactions (02 > 0), is 
the clock model solution with the constant phase differ- 
ence between neighboring layers over the whole sample 
0. The tilt in the j-th layer is £j = 6*{cos(ja), sin(ja)} 
where 9 is constant and the phase difference a is the 
difference in directions of the tilt vectors in neighboring 
layers. The phase difference a is obtained by the mini- 
mization of the free energy Eq. (Q) with respect to a. 

We have analyzed the temperature dependence of solu- 
tions for the following behavior of model parameters with 
decreasing temperature: a\ increases monotonically from 
the negative value to the positive value, piezoelectric pa- 
rameter c p monotonically increases and flexoelectric pa- 
rameter [i monotonically decreases. Let us consider a 
microscopic mechanism of the described parameter de- 
pendence with decreasing temperature. The monotonous 



increase of direct interactions a\ with decreasing tem- 
perature from negative value which is due to interpene- 
trating molecules through nearest layers, to the positive 
value in the region where van der Waals interactions fa- 
vor anticlinic order, is due to increasing smectic order 
p8[ . The flexoelectric coupling {[i) decreases monotoni- 
cally with decreasing temperature, since the smectic or- 
der and hindrance of the rotation becomes less affected by 
the molecules above and below the interacting layer. In 
contrast, with decreasing temperature piezoelectric cou- 
pling (c p ) monotonically increases since higher smectic 
order strengthen the rotation hindrance within the layer. 

For positive ai we can introduce the measure of the 
competition - the competition ratio 45,i/et2. In racemic 
mixtures in the temperature region where the compe- 
tition ratio is smaller than -1, the synclinic nonmodu- 
lated Sm C* phase is stable. In the temperature region 
where the competition ratio is larger than 1, the anti- 
clinic nonmodulated Sm phase exists. In the region 
in between, the phase angle a changes with decreasing 
temperature rapidly from 0° to 180° i.e. modulation 
period changes from infinity to two layers. In slightly 
chiral systems the Sm C* phase appears above the fer- 
roelectric Sm C* phase. With additional increase of the 
enantiomeric excess or increased chirality of the system 
(larger c p ), the ferroelectric Sm C* phase dissapear, leav- 
ing only modulated phase with short pitch (Fig. 1). 

Experimental observations have shown that in racemic 
mixtures only the synclinic Sm C* phase and the an- 
ticlinic Sm phase exist with the first order transi- 
tion between them Q. The narrow temperature range 
of the modulated phase (Fig. 1 - solid line), can be ex- 
perimentally seen as the first order transition. In chi- 
ral samples two additional phases the Sm C* phase and 
the Sm C* phase appear. In chiral systems where piezo- 
electric coupling Cp is important, below the transitions 
temperature to the tilted phase, the modulation period 
is much shorter and increases with decreasing tempera- 
ture. The temperature region where competition ratio 
is between -1 and +1 becomes wider and can be experi- 
mentally observed as a distinguished Sm C* phase. The 
similar behavior is observed in chiral samples of MH- 
POBC f§ and 10OTBBB1M7 || (Fig. 1 - dashed line). 
If piezoelectric c p coupling is still stronger, the ferroelec- 
tric Sm C* phase dissapears, which was observed in some 
systems |l7] (Fig. 1 - dotted line). The present form of 
the model takes into account only quadratic interactions 
and can not account for first order transitions between 
phases with different values of a. 

The symmetry of structures with various a is the same, 
since the symmetry operation common to all structures 
consists from the translation for a layer thickness and 
the rotation for an angle a. If a changes discontinuously, 
the temperature range of different phases can be defined 
with DSC or similar measurements. If transitions are 
continuous, the temperature ranges of various phases are 
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defined by changes of the macroscopic properties and are 
to some extent arbitrary. The phase difference a can be 
used as a parameter which distinguishes between different 
phases. In the Sm C* phase varies from a w 10° to 
70° which means that the modulation period varies from 
40 to 5 layers. In the Sm C* phase is a « i.e. the 
infinite modulation period in racemic mixtures and a few 
hundred layers in pure samples. In the Sm C* FI2 phase is 
a fts 90° or four layers, in the Sm C FI1 phase is a ~ 120° 
or three layers. In the Sm phase is a ss 180°. The 
modulation period is two layers and can also be analyzed 
within a bilayer model where it is treated as a long 
double helix. 

To conclude, we present the phenomenological model 
based on microscopic intralayer and interlayer interac- 
tions. Polarization is induced by piezoelectric and flex- 
oelectric effect and is taken into account explicitly. Al- 
though tilts directly interact only with nearest layers and 
polarizations interact up to the next nearest layers, indi- 
rect achiral interactions which tend to stabilize structures 
with three layer and four layer periodicities extend to the 
fourth neighboring layers. Chiral piezoelectric couplings 
influence also the effective achiral interactions between 
nearest layers. Additionally, chiral interactions due to 
the flexoelectric effect are important to the third near- 
est layers and can be strong. Direct interactions between 
nearest layers correspond to elastic terms in continuous 
models. As a consequence of the piezoelectric and flexo- 
electric coupling, the corresponding Lifshitz parameter in 
continuous models is renormalized. In contrast, variation 
of polarization, which is considered as electrostatic inter- 
layer coupling, has no corresponding terms in continuous 
models. The large variety of phases in antiferroelectric 
liquid crystals is a consequence of a delicate balance be- 
tween various mechanisms. In chiral samples all interac- 
tions change with temperature monotonically. But the 
competition ratio varies nonmonotonically and therefore 
also phase difference a i.e. modulation period changes 
nonmonotonically and the ferroelectric Sm C* phase with 
long modulation period can be stable between two phases 
(Sm C* and Sm C* FI2 ) with much shorter modulation 
periods. 

Within the phenomenological discrete model we stud- 
ied only the simplest possible interlayer interactions and 
the simplest possible clock model solutions. As a future 
problem we left the analysis of the influence of the phase 
difference a on the tilt which could induce first order 
transitions between phases. We also did not consider 
quadrupolar coupling, which could induce quadrupolar 
ordering and therefore allow for distorted clock model 
solutions pg |. 
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FIG. 1. Temperature dependence of the phase difference 
a. Model parameters change with temperature monotonically 
as: oi = (-4.1 K - 3.2 T), c p = x (0.15 K - 0.16 T), 
fj, = (2.12 K + 1.32 T). The other values were / = K; 
bo = 2 K; bi = bo/10; &2 = bo/100. The parameter x gives 
enantiomeric excess and has the following values: x = - 
solid line, x = 0.2 - dashed line and x = 1 - dotted line. 
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